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03 ■ Abstract 

B 

j> | We study the lowest energy E of a relativistic system of iV identical bosons bound 

by pair potentials of the form V(rij) = g(rfj) in three spatial dimensions. In natural 
units h = c = 1 the system has the semirelativistic 'spinless-Salpeter' Hamiltonian 

N N 

i=l j>i=l 

where g is monotone increasing and has convexity g" > 0. We use 'envelope 
theory' to derive formulas for general lower energy bounds and we use a variational 
method to find complementary upper bounds valid for all N > 2. In particular, 
we determine the energy of the iV-body oscillator g(r 2 ) = cr 2 with error less than 
0.15% for all to > 0, N > 2, and c > 0. 



PACS: 03.65. Ge, 03.65.Pm, ll.10.St 



Relativistic N-boson systems 



page 2 



I. Introduction 

We consider a system of N identical bosons interacting by attractive pair 
potentials Vfaj) and obeying the semirelativistic spinless Salpeter equation [1, 2]. 
The Hamiltonian governing the dynamics of the iV-particle problem is given by 

N N 

h = y, ^ 2 +p1 + E v ^ - p ii) ( L1 ) 

i—l jyi=l 

and represents a model system having a relativistically correct expression for the 
kinetic energy and a static pair potential. One of the reasons for considering such a 
model is that the extension to the many-particle case poses no fundamental technical 
problems beyond what are already present in the one-body problem, namely the 
square root in the kinetic energy and the non-locality that the definition [3] of the 
Hamiltonian entails. Our lower bounds use the necessary permutation symmetry 
of the iV-boson problem to effect a 'reduction' to an almost equivalent 2-body 
problem [4,5]. The purpose of the present paper is first to use envelope theory [6-10] 
to extend our specific energy lower bounds for the harmonic oscillator [11] to apply 
to smooth transformations of the oscillator having the general form V(r) = g(r 2 ), 
where g is monotone increasing and of positive convexity ( g" > ) . Secondly, we 
show that the earlier upper energy bounds (via a Gaussian trial function) for the 
oscillator V(r) = cr 2 can be considerably sharpened; this improvement is carried 
over to the larger class of pair potentials. We have already shown this [12] for the 
ultrarelativistic case m = of the pure oscillator. In this paper we shall generalize 
these oscillator results to V(r) = g(r 2 ) and m > 0. For the oscillator V(r) = cr 2 
itself, the new bounds are separated by less than 0.15% for all m > 0, c > 0, and 
N > 2. 

In Section II we recall some fundamental formulas concerning the one-body 
harmonic oscillator with Hamiltonian a/ m 2 + p 2 + r 2 and lowest energy e(m). 
This problem does not have an exact analytical solution but can be easily solved 
numerically to yield e(m) to arbitrary accuracy; this result is necessary for our iV- 
body lower bounds. As distinct from our earlier work [11], in this paper we eschew 
the P-representation and its concomitant scaling subtleties, and base all our lower 
bounds on the function e(m) itself. 

In Section III we turn to the principal topic of this paper, namely potentials 
which are smooth transformations V(r) = g(r 2 ) of the oscillator potential. If g is 
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convex (g" > ), the graph of V(r) lies above 'tangential potentials' V^(r) with 
the general form V^'(r) = a(t) + b(t)r 2 , where t = f 2 is the point of contact with 
the potential V(r) itself. As t > varies, {V^(r)} represents a family of shifted 
oscillators. Envelope theory allows one to construct energy lower bounds based on 
this fundamental geometrical idea. In Section IV we construct variational upper 
bounds by use of a translation-invariant Gaussian trial function. In Section V we 
look at the ultrarelativistic case m — > 0, and in Section VI we apply our general 
results to some examples from the family V(r) = cr q , q > 2. 

II. The one-body oscillator problem 

We consider the one-body problem with Hamiltonian 

H x = A/m 2 + p 2 + r 2 -> e(m), (2.1) 

where, for coupling c = 1, e(m) is the lowest eigenvalue as a function of the mass 
m. In the momentum-space representation, we have an equivalent problem with 
Hamiltonian 

Hi = -A + \/m 2 + r 2 -> e(m). (2.2) 

Since this Schrodinger problem is easy to solve numerically to arbitrary accuracy, 
we shall take the position that e(m) is 'known' and at our disposal. We note that 
in the large-m (nonrelativistic or Schrodinger) limit, we have 

3 

e(m) ~ cnr("z) = m H — . (2.3) 

V2to 

The graph of e(m) — m is shown in Figure 1: e(m) is monotone increasing with 
m; e(m) — m, however, is monotone decreasing , in agreement, for large m, with 
the Feynman-Hellmann theorem for the corresponding nonrelativistic case. 

It remains now to use scaling to generalize these results. This is necessary for 
our later application to the iV-body problem. For the energy of a more general one- 
body problem in which the kinetic-energy term is multiplied by the positive factor 
P, the coupling 7 > is included, and a further parameter A > is allowed for, 
we have, by scaling arguments, 

Hi = P^m 2 + Xp 2 + ~fr 2 -> s(m, (3, 7 A) = (/3 2 7 A) 1/3 e (m (J^J ' J . (2.4) 
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III. Energy lower bound for V(r) = g(r 2 ) by envelope theory 

Our hypothesis is that V(r) = g(r 2 ), where the smooth transformation func- 
tion g is monotone increasing and its convexity is positive or zero. That is to say, 
we shall assume g" > 0. These assumptions imply a relation between V(r) and a 
'tangential potential' V^\r) given explicitly by 

V(r) > V {t \r) = g(t) - tg'(t) + g'(t)r 2 = a(t) + b(t)r 2 , (3.1) 

where t = f 2 is the point of contact between the tangential potential and the 
potential. For each fixed t, the tangential potential has the form a + br 2 of a 
shifted oscillator. This potential inequality induces, in turn, a spectral inequality as 
an immediate consequence of the min-max characterization of the spectrum of the 
Hamiltonian. It is the task of envelope theory [6-10] to generate expressions for this 
spectral inequality. 

The kinetic-energy term in the Hamiltonian H does not have the kinetic energy 
of the center-of-mass removed. Thus the wave function we use must satisfy two 
fundamental symmetries: translation invariance and boson permutation symmetry 
(in the individual-particle coordinates). Jacobi relative coordinates may be defined 
with the aid of an orthogonal matrix B relating the column vectors of the new 
[pi] and old [r*j] coordinates given by [p^ = B[ri\. The first row of B defines a 
center-of-mass variable p\ with every entry l/y/N, the second row defines a pair 
distance p2 = (t*i — r , 2)/v / 2, and the kth row, k > 2, has the first k — 1 entries 
Bki = 1/ yjk{k — 1), the kth. entry B^ = —y/ (k — l)/fc, and the remaining entries 
zero. We define the corresponding momentum variables by [tti] = (B~ ) t [pJ = 
-B[Pj]. Let us suppose that the (unknown) exact normalized boson ground-state 
wave function for the iV-body harmonic-oscillator problem with V(r) = cr 2 is 
^ = ty(p2, P3, ■ ■ ■ , Pn) corresponding to energy E. Boson symmetry is a powerful 
constraint that greatly reduces the complexity of this problem. We immediately 
obtain [11, Eq. (2.3)] the 'reduction' 



E = (#,##) 



NJm 2 +p A N + 



-c\ri — T2\ 



(3.2) 



Since \r\ — r 2 | 2 = 2p2, in terms of the Jacobi relative coordinates this becomes 





( 




E = 




N \ 




V 





m 2 + 



71"! 



N — 1 



N 



-TIN 



+ N(N-l)cp 2 2 



(3.3) 
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The lemma proved in [11] allows us to remove the term in the center-of-mass mo- 
mentum operator tti from inside the square root. Boson permutation symmetry 
furthermore implies [11, Eq. (2.5)] 

(¥,p2*) = (¥,A&¥), (3.4) 

even though the wave function \1/ may not be symmetric in the relative coordinates. 
These results lead to the final reduction 



v v / m 2 + ^_I^ + Ar(A r_ 1)cp 2 r 



¥ . (3.5) 



If we now write r = pn and p = itn, we see that the exact energy E can be 
written in the form E = (^,TC^), in which H is the Hamiltonian for a one-body 
problem given by 

H = P^m 2 + \p 2 +~fcr 2 , (3.6) 

with 

N-l 

p = N, A = ^^, and 7 = iV(iV-l). 

It follows that the exact energy E of the oscillator system is bounded below by £, 
the bottom of the spectrum of the one-body Hamiltonian Ti. 

Thus, for the harmonic oscillator itself, we have from (3.5) and (2.4) 
Theorem 1 

A lower bound to the ground-state energy eigenvalue E of the semirelativistic N- 
body Hamiltonian 

N N 

i=l j>i=l 

is provided by the formula 

E>(lf-,cX) 1/3 e^ra^yy (3.8) 

where 

P = N, A = ^^, 1 = N(N-1). 



This lower bound yields the exact energy in the Schrodinger limit m — > oo. If we 
consider the potential V(r) = g(r 2 ) and use the potential lower bound (3.1), we 
can maximize the resulting lower bound provided by Theorem 1 to obtain 
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Theorem 2 

A lower bound to the ground-state energy eigenvalue E of the semirelativistic N- 
body Hamiltonian 



N TV 

H = J2\J™ 2 +Pl+ E gdn-rjl 2 ), g'>0, g">0, (3.9) 



i=l 



j>i=l 



is provided by the formula 



E > max 

i>0 



m^ + l{g(t)-tg\t)) 



(3.10) 



where 



p = N, A=^^, 7 = iV(iV-l), v = m 



7VW 



1/3 



If we consider the family of pure-power potentials of the form V(r) = cr q , then 
for the harmonic oscillator q = 2, we use Theorem 1; for more general potentials, 
with q > 2, we have V(r) = g(r 2 ) = g(t) = ct q l 2 . Consequently, we must in this 
case make the explicit substitutions: 

a{t) = g(t) - tg'(t) = -c (| - l) t^ 2 and b(t) = g\t) = |^" 2 )/ 2 . (3.11) 



IV. Variational upper bounds 

Improvement over the previous upper energy bounds [11] for the oscillator will 
be obtained in this paper by avoiding the loosening incurred by use of Jensen's 
inequality [3,13]. This goal has already been achieved [12] for the ultrarelativistic 
special case m = of the iV-body harmonic-oscillator problem. We shall now 
extend this to more general problems with attractive potential V(r) and m > 0. 

We use a Gaussian wave function of the form 

®(P2, p 3 , ■ ■ ■ , Pn) =Cexp ^-|E^j ' a>0 ' 

where C is a normalization constant. The factoring property of this function, the 
boson-symmetry reduction leading to (3.5), and the additional fact that $ is also 
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symmetric under exchange of the relative coordinates allows us to write r = p2, 
and p = itn — > 7T2, and finally 

E<p(<f>, Vm 2 + Xp 2 0) + | (</>, ^(v^r)^) , (4.2) 

where 

AT - 1 

P = N, A = ^^, 7 = ^-1), 
and the function 0(r) is given by 

^) = (^) 3/4 e X p(-^). (4.3) 

The kinetic-energy integral may be written in terms of modified Bessel functions 
of the second kind [14,15], which we now discuss. The calculation is best carried out 
in momentum space with the aid of the three-dimensional Fourier transform T . 
We have 

T ( 1 \ 3/4 / P 

— > <j> 

Thus the expectation of the kinetic energy becomes 



«r)^. «*)=(-! exp --). (4.4) 



(K) = (3 ( j>, Vm 2 + Xk 2 $) = - 4n H /0 [ exp ( V 'm 2 + Xk 2 k 2 dk. (4.5) 
V / (any/ 2 J \ a ) 

We may write this integral in the form 

where K v (z) is a modified Bessel function of the second kind [14,15]. 

The potential-energy integral will depend on the choice of V(r). For the family 
V(r) = c sgn((/)r 9 , which we shall study in Section V, the integrals may be expressed 
in terms of the gamma function. Explicitly we have 

With the results in this form we can use the parameter [i as a variational parameter. 
We have therefore established 
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Theorem 3 

For fixed m > 0, q > — 1, c > 0, tV > 2, = N, 7 = iV(iV - 1), and 
A = (N — 1) /A/", i/ie lowest energy E of the N-boson problem for the pair potential 
V(r) = c sgn(q)r q is given by the inequality 

1- 



E < min 

M>0 



exp 



EL 
4 



EL 
4 



c sgn(g)7 r 



3 + <A / A*\/A 



(4.8) 



We have allowed (/ > — 1 here since the upper bound easily accommodates this 
family of potentials in three spatial dimensions. For q < — 1, there is no discrete 
spectrum. In the gravitational case q = — 1 the minimum upper bound exists 
provided the coupling is not too large: specifically we require 



4/3 V A 2 



c 7 2 c N(N-l) 



< 1. 



(4.9) 



This situation is of course well known from the two-particle attractive Coulomb 
problem [16,17]. At present we only have complementary lower bounds for q > 2. 



V. The ultrarelativistic limit 

The ultrarelativistic case m — > may be obtained from Theorems 2 and 3 as 
a special case. The Hamiltonian for this problem is given explicitly by 



TV TV 

H = J2ytf + Yl c ' r * - ^i 9 ' c > 0, ? > 2. 

j>i=i 



(5.1) 



i=l 



For the lower bound we use g 

(£) = c ti/ 2 . The upper bound may either be treated 
separately or taken from (4.8) by means of the limit lhn^o zK\ {£) = 1. The 
bounds we obtain are given by 



C 



£0 
3 



3q si 
*(!+,) < E K C 



2r 



3 + q 



1 

1+9 



(5.2) 



where zq ~ 2.33810741 is the first zero of the Airy function Ai(z), satisfying 
Ai(zo) = 0, and the common factor C is given by 

C = (|) ^ (l + -) (tV(tV - 1))TO 2^T. 
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As m increases from zero, the power-law bounds become closer monotonically with 
m; thus the m = case provides an upper bound to the error for all m. Since we 
have explicit formulas for the bounds in terms of A, we are able to make definite 
statements concerning the percentage separations of the bounds for all A. If we 
take the energy estimate E to be the average of the bounds, the exact energy E 
is determined by E to within 0.15% for the harmonic oscillator q = 2, and to 
3.6% for the cubic 'oscillator' q = 3. In the nonrelativistic limit m — > oo the 
harmonic-oscillator bounds q = 2 coalesce to the exact solution of the Schrodinger 
iV-body problem [11]. 

VI. Examples 

The examples we consider are from the family V(r) = cr q . In order to have 
a lower bound, we restrict the power to q > 2. We revisit the oscillator problem 
because we have considerably improved the upper bound since Ref. [11]. Graphs 
of the lower bounds alone are shown in Figure 2. The percentage separations are 
bounded above by the separations at m = 0, which are there less than 0.15%. With 
the notation Ejsr{m) we have explicitly, for the oscillator q = 2, that the lower and 
upper estimates have numerical values Eg(l) = 35.86383 and Eg (I) = 35.89953, 
respectively. Thus the average of these values determines E 8 (l) in this case with 
error less than 0.05%. As we leave the oscillator and increase q beyond q = 2, 
the bounds become less sharp. For q = | we show the corresponding bounds 
in Figure 3: here the bounds are separated for all m by less than 1.43%. The 
corresponding graphs for the cubic 'oscillator' q = 3 are shown in Figure 4; in this 
case the maximum percentage separation (again for all A > 2 and m > ) is 
3.6%. 

VII. Conclusion 

The necessary permutation symmetry of the states of a system of identical par- 
ticles is a powerful constraint. The approximate 'reduction' of the A-body problem 
to a scaled two-body problem is most striking for systems of bosons, or for systems 
which are compatible with the assumption of permutation symmetry in the spatial 
variables alone [4]. For systems of fermions, the reduction is to a sum over two-body 
energies [5]. For the Schrodinger harmonic-oscillator problem, the boson reduction 
is complete in the sense that the A-body energy is given exactly by the energy of 
a two-body problem. A lower bound by this type of reduction is possible provided 
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either the kinetic-energy term or the potential-energy term has a quadratic form: 
this allows us to replace, for example, the 'mixed' pair {itn, P2} by {^n, Pn} in 
the reduced two-body Hamiltonian H. For the Salpeter problem discussed in this 
paper, a quadratic form is present in the potential term of the oscillator, and the 
lower bound obtained for this base problem is then applicable to other problems 
whose potentials V(r) have the form of smooth convex transformations g(r 2 ) of 
the oscillator. The extension beyond the oscillator is effected by the use of 'enve- 
lope theory'. A reduction is also used for our upper bound, but this reduction is 
allowed for general potentials and for a different reason. The trial function must be 
a translation-invariant boson function; but we have chosen a Gaussian trial func- 
tion which has an additional symmetry, namely, it is also symmetric in the relative 
coordinates. It is this latter symmetry which completes the reduction in the case 
of the upper bound. Because of all these symmetries, what starts out as a complex 
many-body problem, appears in the end, for the purpose of finding energy bounds, 
as a one-body problem. 

Acknowledgements 

Partial financial support of this work under Grant No. GP3438 from the Natural 
Sciences and Engineering Research Council of Canada, and the hospitality of the 
Institute for High Energy Physics of the Austrian Academy of Sciences in Vienna, 
is gratefully acknowledged by one of us [RLH] . 



Relativistic N-boson systems 



page 11 



References 

[1] E. E. Salpeter and H. A. Bethe, Phys. Rev. 84, 1232 (1951). 
[2] E. E. Salpeter, Phys. Rev. 87, 328 (1952). 

[3] E. H. Lieb and M. Loss, Analysis (American Mathematical Society, New 
York, 1996). The definition of the Salpeter kinetic-energy operator is given on 
p. 168; Jensen's inequality is given in Theorem 2.2 on p. 38. 

[4] R. L. Hall and H. R. Post, Proc. Phys. Soc. (Lond.) 90, 381 (1967). 

[5] R. L. Hall, Proc. Phys. Soc. (Lond.) 91, 16 (1967). 

[6] R. L. Hall, J. Math. Phys. 24, 324 (1983). 

[7] R. L. Hall, J. Math. Phys. 25, 2708 (1984). 

[8] R. L. Hall, W. Lucha, and F. F. Schoberl, J. Math. Phys. 42, 5228 (2001). 

[9] R. L. Hall, W. Lucha, and F. F. Schoberl, Int. J. Mod. Phys. A 17, 1931 (2002). 

[10] R. L. Hall, W. Lucha, and F. F. Schoberl, J. Math. Phys. 43, 5913 (2002). 

[11] R. L. Hall, W. Lucha, and F. F. Schoberl, J. Math. Phys. 43, 1237(2002); 
44, 2724 (2003) (Erratum). 

[12] R. L. Hall, W. Lucha, and F. F. Schoberl, Phys. Lett. A 320, 127 (2003). 

[13] W. Feller, An introduction to probability theory and its applications, Volume II 
(John Wiley, New York, 1971). Jensen's inequality is discussed on p. 153. 

[14] M. Abramowitz and I. A. Stegun (eds.), Handbook of Mathematical Func- 
tions with Formulas, Graphs, and Mathematical Tables (Dover, New York, 
1972). Modified Bessel functions of the second kind K n (z) are discussed on 
pp. 374-377. 

[15] Wolfram Research Inc., Mathematica Version 4 (1999). Modified Bessel func- 
tions of the second kind K n (z) are represented in Mathematica by K n (z) = 
BesselK(n,z) . 

[16] I. W. Herbst, Commun. Math. Phys. 53, 285 (1977); 55, 316 (1977) (adden- 
dum). 



Relativistic N-boson systems 



page 12 




Figure 1. The energy function e(m) — m of the one-body problem defined by 
(2.1). 
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Figure 2. The ground-state energy E(m) of the relativistic iV-boson harmonic- 
oscillator problem V(r) = r 2 for N = 2, 3, . . . , 8. The figure shows the lower 
bounds given by Eq. (3.8): the upper bounds are everywhere less than 0.15% above 
these curves and are indistinguishable on the graph. In the Schrodinger limit m — > 
oo the upper and lower bounds coalesce to the exact energies. 
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Figure 3. Upper and lower energy bounds for the ground-state energy E(m) of 
the relativistic iV-boson problem corresponding to V(r) = r 5 / 2 for N = 2, 3, . . . , 8. 
The percentage errors are maximum for m = where they determine the energies 
(for all N) with error less than 1.43%. 
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Figure 4. Upper and lower energy bounds for the ground-state energy E(m) of 
the relativistic iV-boson problem corresponding to V(r) = r 3 for N = 2, 3, . . . , 8. 
The percentage errors are maximum for m = where they determine the energies 
(for all N) with error less than 3.6%. 



